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Phase transitions in a frustrated XY model with zig-zag couplings 
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We study a new generalized version of the square-lattice frustrated XY model where unequal 
ferromagnetic and antiferromagnetic couplings are arranged in a zig-zag pattern. The ratio between 
the couplings p can be used to tune the system, continuously, from the isotropic square-lattice to 
the triangular-lattice frustrated XY model. The model can be physically realized as a Josephson- 
(— I I junction array with two different couplings, in a magnetic field corresponding to half-flux quanta per 

. plaquette. Mean-field approximation, Ginzburg-Landau expansion and finite-size scaling of Monte 

Carlo simulations are used to study the phase diagram and critical behavior. Depending on the value 
of p, two separate transitions or a transition line in the universality class of the XY-Ising model, 
with combined Z2 and U{\) symmetries, takes place. In particular, the phase transitions of the 
standard square-lattice and triangular-lattice frustrated XY models correspond to two different cuts 
through the same transition line. Estimates of the chiral (Z2) critical exponents on this transition 
line deviate significantly from the pure Ising values, consistent with that along the critical line of 
the XY-Ising model. This suggests that a frustrated XY model or Josephson-j unction array with a 
zig-zag coupling modulation can provide a physical realization of the XY-Ising model critical line. 
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I. INTRODUCTION 



^ ■ Thexe has been an increasing interest in frustrated XY models in relation to Josephson-junction arrays in a magnetic 
neldBu. At a particular value of the external field, corresponding to halfrilux quanta per plaquette of the array, the 
ideal system is isomorphic to a frustrated XY model, or Villain's odd models, with ferromagnetic and antiferromagnetic 
bonds satisfying the odd rule, in which every plaquette has an odd number of antiferromagnetic bonds. Frustration 
has the effect of introducing a discrete Z2 symmetry in the ground state with an associated chiral (Ising-like) order 
parameter, in addition to the continuous U(l) symmetry. The interplay between these two order parameters may 
Q\ [ lead to critical behavior which is not present in the unfrustratcd model which is known to have a transition in the 
"p-7 Kosterliz-Thouless (KT) universality class. Ill— ■ n p 
Earlier Monte Carlo simulation results foe, the isotropic square- lattice (SFXY)WO and triangular-lattice (TFXY)ffl 
frustrated XY model, and some recent onesQw, suggest a critical behavior associated with the chiral order parameter 
in agreement with pure Ising exponents while the continuous (XY) degrees of freedom display the main features of the 
KT transition, possibly with a nonuniversal jump. Estimates of the corresponding critical temperatures are always too 
close to be satisfactorily resolved within the errorbars, specially when possible systematic errors due to the assumed 
\ KT scaling forms are taken into account. These results can either be regarded as an indication of a single but decoupled 
transition, where the Ising and XY variables have standard behavior and the same critical point, or else there are two 
separate by close transitions of Ising and KT type. There exist also some appealing arguments which exclude one of 
the two [possibilities, Ising followed by a KT transition for increasing temperature, in the case of a double transition 
scenarioo. Other numerical works, however, which attempt an improved estimate o£ the chiral critical exponents tend 



to conclude that these exponents deviate significantly from the pure Ising valuesollj'EL ptp particular, based on the 
results for the coupled XY-Ising model as an effective Hamiltonian for these systemstJO, it has been argued that, 
in the case of the single transition scenario, both the SFXY and TFXY model display a transition with exponents 
deviating from the pure Ising values. Moreover, the exponents are given bv-the corresponding values along the-critical 
line of this model. Estimates of chiral exponents froauMonte Carlo dataE3 and transfer-matrix calculationsy'tHl are 
consistent with the XY-Ising model universality classEjll3. I— ■ 

A generalized version of the SFXY model has been introduced by Berge et al.EZl where the strength of the an- 
tiferromagnetic bonds can be varied. This introduces a particular anisotropy into the system and leads to clearly 
separated Ising and KT-like transitions for unequal strengthst3~E3 but which appear to merge into a single one for 
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equal strengths, corresponding to the isotropic SFXY model. There is a critical value for the bond strength, 1/3, 
below which the twofold degeneracy disappears. Other generalizations have been introduced for the TFXY model 
that also leads to a critical strength below which the frustration effect is suppressedE3. A common feature in the 
topology of the phase diagram of these generalized versions is that the isotropic model always corresponds to the 
region where chiral and XY ordering can not be clearly resolved. However, so far, the SFXY and TFXY models have 
been treated as separated models. 

In this work, we introduce a new generalized version of the SFXY model where unequal ferromagnetic and antifcr- 
romagnetic couplings are arranged in a zig-zag pattern. The ratio between the couplings p can be used to tune the 
system, continuously, from the isotropic SFXY to the TFXY model, allowing the study of both models within the 
same framework. The model can be physically realized as a Josephson-junction array with two different couplings, 
in a magnetic field corresponding to half-flux quanta per plaquette. We use a mean-field approximation, Ginzburg- 
Landau expansion and finite-size scaling of Monte Carlo simulations to study the phase diagram and critical behavior. 
Depending on the value of p, two separate transitions or a transition line with combined Z2 and U(l) symmetries, 
takes place. Based on an effective Hamiltonian, we show that this transition line is in the universality class of the 
XY-Ising model and the phase transitions of the standard SFXY and TFXY models correspond to two different cuts 
through the same transition line. Estimates of the chiral (Zq) critical exponents are consistent with that along the 
critical line of the coupled XY-Ising model, suggesting a possible physical realization of the XY-Ising model critical 
line in a frustrated XY model or Josephson-junction array with a zig-zag coupling modulation. 

The remainder of the paper is organized as follows. In Sec. II, wc define the model. In Sec. Ill, the ground 
state properties obtained by two different methods are presented. In Sec. IV, a mean field approximation is used to 
obtain the global features of the phase diagram. In Sec. V, the effective Hamiltonian obtained by Guinzburg Landau 
expansions is presented and its relation to coupled XY models and the XY-Ising model is discussed. In Sec. VI, we 
present numerical results of Monte Carlo simulations for the phase diagram and chiral critical exponents obtained 
from finite-size scaling. Finally, Sec. VII is devoted to the conclusions. 



II. THE MODEL 

The generalized version of the frustrated XY model introduced by Berge et alEl can be regarded-as an XY version of 
one of the two frustrated Ising models with periodic interactions first introduced by Andre et al.c3. The other model 
has the important feature that it reduces to the triangular-lattice antiferromagnetic Ising model in one particular 
limit. In analogy to this model, we consider a system of classical XY spins on a square lattice with nearest neighbors 
interactions modulated in a periodic pattern. The Hamiltonian of this zig-zag model is given by 

W = -„ ^ JijSi-Sj, (1) 

<hj> 

where the sum is restricted to the first neighbors and Si is a two-component unit vector. The couplings Jij can have 
two different values, J and J' , distributed periodically in a zig-zag pattern as indicated in Fig. [|. We choose J to be 
ferromagnetic (J > 0) and define J' — —pJ, where p is the coupling ratio. We are interested in the case p > 0, where 
each plaquette has an odd number of antiferromagnet bonds, Villain's odd ruleu, which leads to frustration effects. 
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FIG. 1. Generalized frustrated XY model with zig-zag coupling modulation. Continuous lines correspond to coupling Jy = J 
and double lines to = J'. The sites A and B denote two sublattices where spins do not interact. 
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When p = 1, the model reduces to the SFXY model while in the limit p — > +00 it is topologically equivalent to 
the TFXY model. The latter limit can be easily established after performing the "gauge" transformation Si — > e^, 
Jij — > JijEiej, where e, = 1 and —1 on the sublattices A and B of Fig. |], respectively, resulting in antiferromagnetic 
J and ferromagnetic J' couplings. When p — > +00, each pair of spins connected by a J' bond become locked and 
may be replaced by an effective spin, leading to an antiferromagnetic XY model with the same coordination number 
as the triangular lattice. The model is then well suited for the study of the universality classes of both SFXY and 
TFXY models. When p — 0, Eq. ([j]) reduces to a ferromagnetic XY model on a hexagonal lattice which undergoes a 
KT transition. 



III. GROUND STATE 

In another generalization of the frustrated XY modeJll, it has been shown that the lowest energy state can be 
constructed by building up the configuration of the infinite lattice from the ground state configuration of a single 
plaqucttc. In our case, the same procedure can be used if we allow for rotations and reflections of the "one plaqucttc 
ground state configuration" which also assures that the true ground state is obtained. No assumption on the periodicity 
of the ground state is made. The plaquette configuration is indicated in Fig. [2]a for p > 1/3 and it is the same as used 
111 Ref. [P?]. The spin configuration is collinear for p < 1/3 and a canted one for p > 1/3. For the canted configuration 
one can define a chiral variable for each plaquette P 



^ ' JijSi x Sj, (2) 
a ° <ij>eP 

where X)<y>eP * s a direct summation around the plaquette and 0~ o is a normalization factor given by 



<ij>£P 



3P+1 3p-l 

For p > 1/3, the ground state of the infinite system constructed by the above procedure consists in a helical spin 
ordering which is incommensurate with the underlying square lattice, except when p — 1 and p = 00, corresponding to 
the SFXY and TFXY models. The pitch A of the helical configuration, can be obtained from half the phase difference 
within the same sublattice in the x direction and is given by 



A = 2cos- 1 (ij^i^) = arccos( ^— 9 - ) (4) 



2 V P V 2 P 

In Fig. |b we show the resulting ground state configuration. The ground state is double degenerate, corresponding 
to an antiferromagnetic arrangement of plaqucttc chiralities a — ±1. 




FIG. 2. Ground state for p > 1/3 consisting in a helical configuration of spins, a) ground state configuration of a single 
plaquette and b) spin configuration for the infinite lattice. 
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For p < 1/3, where the single plaquette configuration is collinear, the ground state is a ferromagnetic configuration 
of spins. 

As an alternative to the above method, the ground state can also be obtained by a direct minimization of the Fourier- 
transform interaction matrix, —Jq, q >. In the present case, we note that there are two non-interacting sublattices, 
corresponding to the sites A and B in Fig. [j], where the Fourier transform can be easily carried out. The interaction 
matrix J q , where k, I denote the sublattices A and B, can be written as 

ju _ , / e iq * - pe~ iq * + e iq y + e~ 1 ^ 

i ~ \ -pe iq * + er iq * + e tq y + e~ iq y 

The eigenvalues are given by, X q — ±V q , where 



V q = J \J i 1 + Pf + 4(-pcos 2 (g x ) +cos 2 (q y ) + (1 - p) cos(q x ) cos(q y )), (6) 

and the dominant eigenvalue X q — +V q reaches a maximum at 

(Qx = 0,q v = 0) for p < i , , 

(q x = ± arccos [(1 - p) / (2p)] , q y = 0) for p > | 1 ' 

From Eq. (|7|), the wavevector characterizing the ground state for p > 1/3 is, in general, incommensurate with the 
lattice periodicity in the x direction, except for p — 1 and p — oo corresponding to the SFXY and TFXY models. The 
eigenvector associated to the largest eigenvalue is a possible realization of ground state provided the corresponding 
spin configuration satisfy the unit vector condition \Si\ = 1. In the present case, they do satisfy this condition and 
correspond to the same configuration as found from the single plaquette method described above. 




IV. MEAN-FIELD PHASE DIAGRAM 



The general form of the phase diagram can be obtained by a mean field (MF) analysis. Although, at finite 
temperatures, this analysis neglects the role of fluctuations it gives nevertheless a good qualitative picture of the 
phase diagram that can also be greatly improved by perturbative or variational techniques. The details are described 
in the Appendix. 
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FIG. 3. Mean field phase diagram. F indicates the ferromagnetic phase, P the paramagnetic phase and H the helical (chiral) 
phase. 

Fig. H shows the phase diagram obtained by the mean field approximation. For p < 1/3, the system undergoes 
a transition from a paramagnetic to a ferromagnetic phase along the transition line XL. This transition is in the 
KT universality class since there is only a single critical mode q = (0, 0) and no additional symmetry in the ground 
state besides the continuous U(l) symmetry. For p > 1/3, there is a paramagnetic phase at high temperatures 
and a helical phase at low temperatures which is incommensurate with the lattice periodicity except for p = 1 and 
p + oo where the model reduces to the SFXY and TFXY models, respectively. The helical phase has an additional 
discrete Zi symmetry associated with the antiferromagnetic arrangement of plaquette chiralities ap in the ground 
state. The mean field analysis gives a single transition for p > 1/3. Therefore, the whole line LT for p > 1, including 
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the TFXY limit, is expected to have the same type of behavior as the SFXY model at p = 1. The nature of this 
transition however can not be studied at mean field level and other methods are required, as will be presented in See. 
V and VI. Note that, in contrast to the generalized version of the frustrated XY model considered by Berge et aZ.tZI, 
where a clear separation into two transitions for p ^= 1 is already found at the mean field leveO, the zig-zag model 
displays two transitions only for p < 1/3 within the same kind of approximation. This suggests that the separation 
of these transitions is not simply a result of the induced anisotropy for p ^ 1 but should be related to the nature of 
the coupling between chiral (Ising) and XY degrees of freedom, in agreement with arguments based on an effective 
coupled XY-Ising model Hamiltoniaro. As will be shown in Sec. IV, for the zig-zag model the form of this coupling 
is unchanged for p ~ 1 and p > 1, suggesting that a clear separation is not expected. 

The transition line CL separating the ferromagnetic from helical phase can be regarded as a commensurate- 
incommensurate transition which joins the other transition lines, XL and LT, at a Lifshitz point L at T / 0. In 
mean field, this transition line is given by p = 1/3 corresponding to the stability boundary between the two modes in 
Eq. (Q). Although, there are interesting .questions regarding the precise location of the Lifshitz point and the nature 
of the phase transition along this lineSTEa, these will not be the subject of a detailed study in this work. 



V. EFFECTIVE HAMILTONIAN 



The universality class of phase transitions can be considered on the basis of an effective Hamiltonian obtained by 
Ginzburg-Landau expansions. Invoking the universality hypothesis, one expects that models with the same effective 
Hamiltonian differing only by irrelevant terms are in the same universality class. In this section, we discuss the critical 
behavior in the region p > 1/3 where the ground state is double degenerated by deriving the corresponding effective 
Hamiltonian. An effective Hamiltonian can be obtained from the free energy functional, describing fluctuations around 
the MF solution discussed in Sec. IV, via a Hubbard- Stratonovich transformation in a standard wayEJEj. One replaces 
Eq. © by 

^ = \H K ^-^-H w ^ (8) 

where ij are unconstrained spins weighted by W(x) ~ x 2 /4 — x 4 /64 + 0(x 6 ) and Kuj = Jij/kT. In the present case, 
we can separate the lattice spins into two non-interacting sublattices, corresponding to sites A and B in Fig. [l| The 
interaction matrix J l q ' k is then given by Eq. (|^) and the corresponding eigenvalues by Eq. (JsJ) . For p > 1/3, there are 

two degenerated modes (fig- and 4>q+, that maximizes the dominant eigenvalue A = +Vq, corresponding to the wave 
vectors in Eq. (0). Retaining these modes only and introducing the real two-component fields 4>i — \{$q- + 4>Q+) 
and (f>2 = 5«(0q- — 0Q+) i one can expand Eq. (||) to quartic order in 0i,2 leading, in the continuum limit, to a free 
energy density of the form 

+ +»((<?!• $if - $4), (9) 

where r = kT/\Q — 1/2, e = j^r-^p-^Q , f — ^r-^i^Q , an d u, v > 0. For p < 1/3, there is only one critical 

mode, (q x ,Qy) = (0,0), and the resulting Ginzburg-Landau expansion has a single two-component fluctuating field 
which is known to lie in the KT universality class. Apart from the space anisotropy, e ^ f when p ^ 1, that can be 
eliminated by rescaling the x and w-spaeej-directions appropriately, the free energy (^|) has the same form as those 
obtained for the SFXY and TFXYE30E2H13 in terms of complex scalar fields ip% — \ip\e iei . In particular, since the 
present model incorporates both the SFXY and TFXY as special cases, it clearly demonstrates that the SFXY and 
TFXY are described by the same Ginzburg-Landau free energy up to quartic order, in agreement with the arguments 
ofRef. p. 

As usual, in two dimensions, fluctuations in the magnitude of the order parameter are assumed to be irrelevant. 
We can then approximate these magnitudes by their corresponding mean-field values {tpi^l = ipo = —r {2u — v) and 
consider only fluctuations of the phase 9i in Eq. (|^) , leading to an effective lattice Hamiltonian in the form of two 
coupled XY models 

[3H = -J2 [ r icos(0i,i - 9 ld ) +r 2 cos(0 2 ,; - 9» d )] - /i^cos2(0 M - (10) 

<ij> i 
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where Ti = T2 = \ip \-\fej and thej-spatial anisotropy has been removed by rescaling x — > x^J ejf, y — > y. In a 
renormalization study of this modeltJ, the T\ = T2 subspace is only preserved under renormalization if they are 
initially equal. For Ti ^ T2, a double transition is found with an Ising followed by a KT transition as temperature is 
increased. It is also found that the coupling term ft, is a relevant variable locking the phase difference into 621 = On 
where r = 0, 1. This leads, in the h — > 00, to an effective Hamiltonian in the form of coupled XY and Ising modelai.ll 

f3H=~Y J [iAefi + B eJ jtT i( Tj)cm{9i - 9j) + CeffViaj], (11) 
<ij> 

where A e ff, B e ff and C e // are effective couplings which depend on the initial values of I\2, h and other couplings 
generated by the renormalization procedure, and o~i = 2^—1 = ±1 is an Ising-like variable. The condition A e ft = B e ff 
is preserved if Ti = I^ in Eq. ([h]) as is the case for the zig-zag model, even though thisjjnodel is anisotropic for p ^ 1. 
This should be contrasted to the generalized SFXY model considered by Berge et alLJ where the Ginzburg-Landau 
free energy has the same form as in Eq. (^J) but with a spatial anisotropy, in the x and y directions, of different 
magnitudes for the cf>i and $2 fields which can not be removed by simple rescalingEl. This leads to coupled XY models 
with Ti ^ T2 in Eq. (|l0|) and consequently should be described by an XY-Ising model with A ^ B in Eq. ( |TT| ) which 
undergoes two seoarate transitions, an Ising followed by a KT transition for increasing temperatures, in agreement 
with simulationsEZI. 

The phase diagpam of the XY-Ising model of Eq. ( |TT| ) for A = B consists of three branches which meet at a 
multicritical pointlij. One of the branches corresponds to single transitions with simultaneous loss of XY and Ising 
order, and the other two to separate KT and Ising transitions. The line of single transitions eventually becomes first 
order further away from the branch point. Our model corresponds to a particular path through the phase diagram 
of the XY-Ising model and the single or double character of the transition depends, fm the relative position to the 
multicritical point. Since there are already indications from numerical simulationsall3'E3 that both SFXY and TFXY 
limits are in the single transition region, we expect that the whole transition line LT for p > pL in Fig. [}] should 
correspond to this critical line. Numerical estimates of critical exponents associated with the Z 2 order parameter for 
the XY-Ising model deviate significantly from the pure Ising values along the critical lineEMLa and will be used in Sec. 
VI to identify which particular path through the phase diagram is realized for the zig-zag model. 

VI. MONTE CARLO SIMULATIONS 

Due to the presence of an incommensurate phase, the standard periodic boundary conditions are not appropriate for 
the zig-zag model since they cause an additional frustration in the system. Therefore we use a self-consistent bmmAary 
condition that allow the system to adapt the boundary condition to the pitch of the helical configuraticmEjEZI. In 
addition, this boundary condition also improves the determination of the spin stiffness. A similar methods has also 
been used recently for the SFXY. 

A. Phase diagram 

To determine the global phase diagram we used simulations of a 36 x 36 system for various values of p. For each 
value, two separate simulations, one starting from the ground state and the other from the high temperature phase, 
were used to estimate the critical temperature. The transition temperature, Ti, associated with the chiral order 
parameter, was obtained from the peak in the chiral staggered susceptibility, with the chiral order parameter defined 
by Eq. (^|). An estimate of the KT transition temperature, Tkt, was obtained from the expected universal value of 
the spin stiffness 7, ^{!Tkt)I^Tkt — 2/7r , at the transition. Since for p ^ 1 the model is anisotropic, 7 was obtained 
as 7 = ytjxxjyy , where j xx and "f yy are the x and y components of the stiffness 7y. This is the quantity that should 
be universal at the KT transition. The anisotropy of the renormalized Gaussian model, at the critical point, can 
be removed by rescaling the coordinate axes, as in the derivation of Eq. ( |l0|), leading to the geometric mean as the 
effective stiffness. The same averaging procedure has also been used in Ref. |20] to treat the anisotropic model of Ref. 

The phase diagram obtained by Monte Carlo simulations is shown in Fig. W. The estimates of Tj and Tkt agree 
within the errorbars for p larger than a critical value p^, which we take as an estimate of the Lifshitz point. This phase 
diagram is similar to the MF result of Sec. Ill but the Lifshitz point is located at (pl — 0.6, Tl — 0.35) and should 
be compared with the MF result, \p = \,T = ?). The transition line XL has the main features of a KT transition 
with a jump in the spin stiffness consistent with the universal value 2/7T and a nondivergent specific heat. 
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FIG. 4. Phase diagram obtained from Monte Carlo simulations. For p > pi, critical-temperature estimates Ti and Tkt agree 
within errorbars and only Ti is indicated. 

The transition line CL is characterized by a divergent chiral susceptibility and an apparent continuous vanishing of 
j xx while 7j,j, r£tnains finite as shown in Fig. [i| An analysis similar to the one used in Ref. for the model studied 
by Berge et al.EZI can also be applied to the zig-zag model and shows that "f xx is inversely proportional to the chiral 
susceptibility and should therefore decrease continuously at the transitio^i-whcn the susceptibility diverges. Similar 
behavior has been found in a generalized model for the triangular latticeoull. 



2.0 




FIG. 5. Temperature dependence of spin stiffness and chiral susceptibility through the CL transition line of the phase diagram 
in Fig. [| at p = 0.4. The data points for ^ xx are scaled by 10 and for x° scaled by 1/10. 



B. Critical exponents 



There have been eetjently many attempts to obtain improved estimates of the critical exponents for the fully 
frustrated XY modeOat3. For the continuous symmetry, the available scaling forms requires the simultaneous fit 
of two or more parameters and an assumption of KT behavior. This may lead to systematic errors in the location of 
the KT transition temperature. For the chiral ( Ising-like) order parameter there exist scaling analysis which do not 
require a precise knowledge of the bulk T c and can provide an estimate of the critical exponents with a one-parameter 
fit. As simple estimates of T/ and Tkt already agree within errorbars along the transition line for p > p^, as indicated 
in Fig. U, attempting to locate the transition line using, separately, KT scaling forms for the U(l) symmetry and pure 
Ising critical behavior for the chiral variables will inevitably lead to estimates of critical points which are difficult to 
resolve on purely numerical grounds due to errorbars. However, if the critical behavior along this line is in fact in the 
same universality class as the XY-Ising model as suggested by the analysis-of Sec. V, then in order to verify the single 
nature of the transition, it is sufficient to study the Z 2 degrees of frccdomlla. If the critical exponents are inconsistent 
with pure Ising values, the transition cannot correspond to the Ising branch of a double transition or to a single but 
decoupled transition. Moreover, the value of the critical exponent can be used to verify if indeed the critical behavior 
corresponds to the critical line of the XY-Ising model. In order to estimate the chiral critical exponents independently 
of T c , we use the same method, based on the finite-size scaling of free energy barriers, which has been applied to the 
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SFXY and TFXY modelsEl. 

In order to obtain good statistics, we consider only systems of size 8 to 36 x 36, with typically 6 - 12 x 10 6 Monte 
Carlo steps. The simulations were performed near the effective (finite size) critical temperature found in the previous 
Section. The histogram method is then used to extrapolate the needed quantity for different temperatures in the 
vicinity of the critical temperature. We follow the same method used in Refs. [l5p| for the SFXY and TFXY models. 
The thermodynamic critical temperature T c can be determined by the crossing of the free-energy barriers AF(T, L), 
obtained from the chirality histogram N(a) as AF = Am(T, L) — A m (T, L), where Am is the maximum and A m one of 
the two minima in A(a) = — log N(a). At the critical point AF is scale invariant but sufficiently close to T c , AF can 
be expanded to linear order in tL x l v as AF = a + btL 1 ^, where t = (T — T c ) /T c . In this scaling regime, the exponent 
vj can be extracted from the finite-size behavior of the temperature derivative = bL x l v as a one-parameter fit 
in a log-log plot, without requiring a precise (simultaneous) determination of T c . The exponent 1(5 jv is extracted 
from the scaling behavior of <7 m j„ ~ L~®l h ', corresponding to the minimum A m (a ln i n ), which only holds at the critical 
point and thus is more affected but the estimate of T c . 
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FIG. 6. Finite-size scaling of the free energy barrier AF for p — 1.5. 

We have studied two different values of p in detail, p = 0.7 and p = 1.5, which are located between the SFXY model 
limit and the Lifshitz point and between the SFXY model and the TFXY model limit, respectively. For p = 1.5 we 
observed crossing of AF for L > 18 as shown in Fig. ^. Corrections to scaling are clearly seen for 6 < L < 12. 
These sizes were not used for the estimates o£-Gritical exponents. Note that, this free energy barriers suffer less from 
corrections to scaling than Binder's cumulanttHl, Ul = 1— < a\ > /Z < a\ > 2 , which is also expected to cross at a 
unique point. This is shown in Fig. ^ where a sign of unique crossing is only observed at the largest system sizes. The 
latter behavior has been used by Olssorfl in relation to the SFXY model to suggest that there are in fact two separate 
transitions and the estimates of v are still dominated by small system sizes. The method we are using, however, 
indicates clearly a single crossing point suggesting a reliable estimate of vi. Fig. || shows the size dependence of the 
slope ^|t c from where l/v = 1.25(1) can be estimated and Fig. ^| shows the behavior of cr m j n which gives the 
estimate 2/3/ 1/ = 0.29(2). The critical temperature is obtained from the value of T at the crossing point in Fig. |?] 
and gives T c = 0.4935(5). The same analysis has been done for p = 0.7 giving T c = 0.408(2), l/v = 1.28(2) and 
2(3 /v — 0.32(4). These estimates deviate significantly from the pure Ising values, l/v = 1 and 2(3/v = 1/4, and suggest 
a single trapsition scenario. Moreover, these values are consistent with those found for the XY-Ising model along the 
critical lineEjEJ. This is the same behavior found for the SFXY and TFXY models using the same methodst3. 
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FIG. 7. Finite-size scaling of Binder's cumulant Ul for p — 1.5. 
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FIG. 9. Finite-size scaling of a m in for p = 1.5. 
VII. CONCLUSION 

We have introduced a new generalized version of the square-lattice frustrated XY model where unequal ferromagnetic 
and antiferromagnetic couplings are arranged in a zig-zag pattern. One of the main features of the model is that the 
ratio between the couplings p can be used to tune the system through different phase transitions and in one particular 
limit it is equivalent to the isotropic triangular-lattice frustrated XY model. The model can be physically realized as a 
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Josephson-junction array with two different couplings and in a magnetic field corresponding to a half-flux quanta per 
plaquette. We used a mean-field approximation, Ginzburg-Landau expansion and finite-size scaling of Monte Carlo 
simulations to study the phase diagram and critical behavior. Mean-field approximation gives a phase diagram which 
qualitatively agrees with the one obtained by Monte Carlo simulations. Depending on the value of p, two separate 
transitions or a transition line with combined Z2 and U (1) symmetries, takes place. Based on an effective Hamiltonian, 
we showed that this transition line is in the universality class of the XY-Ising model and the phase transitions of the 
standard SFXY and TFXY models correspond to two different cuts through the same transition line. Estimates of 
the chiral (Z2) critical exponents from a finite-size analysis of Monte Carlo data were found to be consistent with 
previous estimates for the SFXY and TFXY models using the same methods. They also agree with the corresponding 
values along the critical line of the coupled XY-lsing model suggesting a possible physical realization of the XY-lsing 
model critical line in a frustrated XY model or Josephson-junction array with a zig-zag coupling modulation. 



APPENDIX: 



The mean field equations for the zig-zag model can be derived by an analysis similar to the one used in Ref. [M 
The corresponding MF equations are 



Mi = S.. 



MF 



R(PH 



Hi 



where H = J2j JijMj is the mean field, R(x) — I\(x) / Iq(x) and (3 = 1/kT 



for 



0, which reduces to 



To find the MF phase diagram we expand (Al) about the transition temperature T c using R(x) 



M i=\^FzZ J ^' 



(Al) 



0(x) 



(A2) 



It appears that one needs to make an assumption on the form of the solution Mj in order to find T C MF . However, if we 
note the similarity of Eq. ( A2) and the zero temperature limit of Eq. ( Al), we can identify the transition temperature 
as 



rpMF 



GS 



(A3) 



provided the local field H L = Jy Mj is independent of the position. Although this property is not expected to hold 
in general, it is satisfied exactly in the ground state found in Sec. III. We then obtain 



> 



T. 



MF 



rpMF _ 



(-P + 3J/2 

\/(l + p) 3 /4p 

If, in addition, we assume that Hi remains independent of i at any temperature < T < we obtain 

R(/3H) 1 



(A4) 



H 



H GS- 



(A5) 



This equation, together with Eq. (Al), shows that the structure of the local configuration around a plaquette and 
the pitch of the helical configuration is independent of the temperature in this approximation. 

For p — > +00 we expect to retrieve the mean field solution of th e TFXY model. However Eq. (A4) leads to a 



diverging value of T c as p — ► +00. As can be seen from Eq. (Al), the temperature is scaled by the magnitude of 



the mean field vector ( J2j JijMj ) which diverges when — > +00. This is an artifact of the mean field approximation 

and other methods, such as perturbative or variational approximationB , can remove this divergence. In fact, the 
phase diagram obtained by Monte Carlo simulations in Sec. V leads to a transition temperature that saturates, for 
p — * +00, to a value consistent with the transition temperature of the TFXY model. 
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